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SUMMARY 


PROBLEM 

Develop  a  mathematical  model  to  predict  the  scattering  of  a  plane  acoustic  wave 
from  an  elastic  spherical  shell  with  different  fluid  media  inside  and  out. 

RESULTS 

A  mathematical  model  of  the  spherical  she!!  was  developed  from  the  more  conveni¬ 
ent  standpoint  of  elastic  theory  rather  than  from  the  accnstic  point  of  view.  The  resulting 
equations  were  solved  in  terms  of  known  functions.  Matrix  methods  are  used  throughout. 
A  numerical  method  is  presented  that  takes  advantage  of  fjme  peculiarities  of  the  resulting 
matrix  and  is  appropriate  to  computer  implementation.  Som?  results  of  a  computer  pro¬ 
gram  for  the  spherical  shell  are  given.  A  few  limiting  cases  of  the  exact  solutions  are  com¬ 
pared  with  previously  reported  results. 
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INTRODUCTION 

In  this  paper,  a  mathematical  model  is  developed  to  predict  the  scattering  of  a  plane 
acoustic  wave  from  an  elastic  spherical  shell  with  different  fluid  media  inside  and  out.  The 
model  is  developed  from  the  more  convenient  standpoint  of  elastic  theory  rather  than  from 
the  acoustic  point  of  view.  The  resulting  equations  are  solved  in  terms  of  known  functions. 
Matrix  methods  are  used  throughout.  A  numerical  method  is  presented  that  is  appropriate 
to  computer  implementation  and  takes  advantage  of  some  peculiarities  of  the  resulting 
matrix.  Some  results  of  ?  computer  program  for  a  spherical  shell  are  given.  A  few  limiting 
cases  of  the  exact  solutions  are  compared  with  previously  reported  results. 

The  paper  consists  of  two  sections  in  addition  to  the  Introduction.  The  first,  A 
Summary  of  Elastodynamics,  summarizes  the  essential  elements  of  elastodynamics  needed 
to  provide  a  unified  treatment  of  the  scattering  problems  considered  in  this  paper.  The 
second.  Scattering  from  a  Spherical  Shell,  develops  the  desired  mathematical  model. 


SUMMARY  OF  LINEAR  ELASTODYNAMICS 

Problems  in  elastodynamics  are  similar  to  boundary-initial  value  problems  in  other 
areas  of  mathematical  physics  and  arc  characterized  by  a  fundamental  system  of  field  equa¬ 
tions  with  corresponding  initial  and  boundary  conditions.  Given  a  prescribed  system  of 
forces  and  displacements  over  the  surface  of  an  elastic  body,  we  may  predict  the  deforma¬ 
tion  of  the  body  using  a  solution  of  the  appropriate  elastic  problem.  Solutions  to  the  elastic 
problems  considered  in  this  paper  are  constructed  using  scalar  and  vector  potentials.  The 
basic  laws  of  infinitesimal  linear  elasticity  and  the  relevant  elements  of  potential  theory  are 
presented  here  for  reference. 

EQUATIONS  OF  MOTION 

The  law  of  motion  in  linear  elasticity,  usually  referred  to  as  Cauchy’s  law,  is  given  by 

divT  =  pu  ,  ( 1 ) 

where  T  is  the  stress  tensor,  u  is  the  displacement  vector,  and  body  forces  have  been  ne¬ 
glected.  In  this  paper  we  will  assume  infinitesimal  displacement  theory  is  valid  so  that  the 
strain  tensor  E  is  related  to  the  displacement  by 

E  =  \  <yu  +  VuT) .  (2) 

In  addition,  we  will  restrict  our  attention  to  homogeneous  and  isotropic  materials  where  the 
constitutive  relation  is  given  by 


I 


I 


) 


T  =  2pE  +  A(trE)  I  , 


(3) 


where  X  and  p  are  the  Lame  constants  and  I  the  identity  tensor.  Using  a  few  standard 
identities,  equations  1  to  3  may  be  manipulated  to  obtain  Navier’s  equation  of  motion: 

/iV2u  +  (X  +  ju)V(V*u)  =  p u  .  (4) 

This  equation  may  also  be  written 

Cl^V<V  *u)  -  Cy2  curl  curl  u  =  u  ,  (5) 

where  the  longitudinal  and  transverse  elastic  wave  velocities  are  given,  respectively,  by 

(6) 

INITIAL  AND  BOUNDARY  DATA 

In  the  mixed  problem  of  elastodvnamscs  we  must  specify  the  following  initial  and 
boundary  data  to  insure  a  unique  solution: 


initial  conditions 

u(*,0)  =  uQ  and  u(*,0)  =  j^onB 

(7) 

Displacement  boundary  condition 

u  =  u  on  Sj  X  (o,t) 

(8) 

Traction  condition 

t  =  jt  on  St  X  (o.t) 

(9) 

where  uQ,  v0.  u,  t  are  specified  data.  and_B  is  an  open  set  in  Euclidean  space  £  which  is  the 
interior  of  a  closed  and  bounded  region  B  (the  closure  of  B)  whose  boundary  3B  is  the 
union  of  surfaces  dB  =  Sj  U  S2. 

POTENTIAL  THEORY 

Next,  we  will  briefly  outline  the  basic  elements  of  potential  theory. 

Theorem  A 

Let  4>  be  a  scalar  field  that  is  continuous  on  B  and  class  on  B.  Then  there  exists 
a  scalar  field  ysuch  that 

<t>  =  V~y  .  (10) 

Proof'  Sec  Kellogg  (reference  1 ) 

The  above  theorem  is  easily  extended  to  vector  fields  by  applying  the  result  to 
each  component. 


Theorem  B  -  Helmholtz  Theorem  (reference  2) 

Let  u  be  a  vector  Held  which  is  continuous  on  B  and  class  on  B.  Then  there 
exists  a  scalar  field  0  and  a  vector  field  0,  both  of  class  C^,  su.h  that 

u  =  V0  +  v  X  0  and  div  0  =  0  . 

Proof:  From  Theorem  A,  there  exists  a  vector  field  v  on  B  such  that 
u  =  V-v  . 

Using  the  standard  identity  curl  curl  v  =  Vdiv  v  -  }  we  fini» 

u  =  vdiv  v  -  curl  curl  v  . 

Therefore,  if  we  let 

4>  =  div  v  and  0  =  -  curl  v  , 
the  desired  conclusion  follows. 

Theorem  C  -  Lame'  Solution  (reference  2) 

Let  u  =  V0  +  V  X  0,  where  0  and  0  are  class  fields  on  B  that  satisfy 

V~  <t> - ^0  =  0 


V20--L0  =  O.  ( 

cT-  ~ 

Then,  u  is  an  clastic  motion  which  satisfies  Navier’s  equation  (equation  5). 

Proof: 

Cp"  Vdiv  u  -  Op2  curl  curl  u  -  U  =V  |cp2v“0-0j  +curl 

jj-  Cj“  curl  curl  0  -  0J  =  V  j^Cp2 V  2  0  -  0-J  +  curl  Jc-p2  V“  0  -  =  0 

where  curKcurl  curl  0)  =  -  curl(v2  0)  by  the  previous  identity. 

The  Lame  solution  can  be  shown  to  be  complete,  that  is,  any  solution  to  Navier’s 
equation  can  be  written  as  a  scalar  and  vector  potential.  The  completeness  theorem  and 
proof  arc  discussed  in  Sternberg  (reference  2). 

Notice  that  the  potentials  0  and  0  are  not  unique.  For  example,  let 

0'  =  0  + constant  0'  =  0  +  vy  . 

The  displacement  u  is  unchanged  using  either  the  primed  or  unprimed  potentials. 


FLUID  MEDIA 

In  this  paper  a  fluid  Will  be  treated  as  an  elastic  material  which  supports  no  trans¬ 
verse  waves,  Cy  =  p  =  0.  In  this  case  Navier’s  equation  becomes 


C2V(V*u)  =  u  c2  =  p-  (16) 

Using  the  Lame  solution  (Theorem  C,  \p  -  0) 

u  =  V0,  whereby  equation  14,  (17) 

V20-4,0  =  O.  (18) 

C2 

The  stress-displacement  relation  may  be  written 

T  =  pC2(divu)I  =  pip  I  .  (19) 


In  many  cases  the  time-dependent  solution  can  be  written  as  a  superposition  of 
sinusoidal  solutions,  for  example  using  Fourier  Transforms.  If  we  assume  the  solution  is 
smooth  enough  that  we  can  interchange  the  differentiation  and  transform  operations,  we 
may  construct  the  time-dependent  solution  frequency  by  frequency.  Throughout  this  paper 
we  will  assume  the  functions  to  be  smooth  and  their  time  dependence  to  be  c' ',on.  For 
convenience  we  will  summarize  the  basic  equations  discussed  in  this  chapter  removing  the 


time  dependence. 

SUMMARY  OF  SOLID  MATERIAL  EQUATIONS 

Na'.^r’s  equation 

— ^V2u - ^  curl  curl  u  =  -u  (20) 

kL~  ~  kT“  ~  ~ 

Lame  solution 

u  =  Vtf>  +  V  X  \p  (21) 

V20  +  kL2^  =  O,  kL2  (22) 

CL~ 

V-iiz+ky2!/'  =  0,  ky-  ~  03  2  (221 

~  CT 

Stress  displacement 

T  =  p [vu  +  VuTl  +  X(div  u)  I  (24) 
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SUMMARY  OF  FLUID  MATERIAL  EQUATIONS 
Navier’s  equation 


V(V’u)  =  -k2u 

k2  = 

C2 

(25) 

Lamd  solution 

u  = 

(26) 

V2  <fi  +  k2  0  =  0 

(27) 

Stress  displacement 

T  =  -p  co2  <j>  I 

(28) 

SCATTERING  FROM  A  SPHERICAL  SHELL 

INTRODUCTION 

A  sphere  is  one  of  the  few  finite  geometrical  shapes  in  which  the  equations  of  linear 
elasticity  can  be  solved  analytically  in  terms  of  known  functions.  Such  a  solution  is  useful 
not  only  for  direct  sonar  mathematical  modeling  applications,  such  as  spherical  baffles, 
domes  and  acoustic  lenses;  but  serves  as  a  useful  cheek  on  various  approximate  techniques. 
For  example,  a  comparison  of  resonant  frequencies  is  made  with  various  thin-shell  theories 
in  reference  3.  Because  the  solution  is  exact,  experimental  verifications  may  be  obtained 
for  a  large  range  of  parameters,  the  limits  being  upon  the  experimental  apparatus. 

In  this  section,  a  solution  suitable  for  digital  computer  implementation  is  developed 
for  plane  wave  scattering  from  an  elastic  spherical  shell  with  different  fluid  media  inside  and 
out.  The  problem  for  identical  fluids  inside  and  outside  has  been  done  by  Goodman  and 
Stern  (reference  4).  In  addition,  numerical  implementation  details  and  a  few  limiting  cases 
are  discussed. 

THEORETICAL  DEVELOPMENT 

Consider  a  homogeneous  plane  wave  impinging  upon  an  elastic  spherical  shell  of 
outer  radius  a,  inner  radius  b  and  thickness  5.  Let  the  incident  wave  come  in  along  the 
polar  axis  (with  no  loss  in  generality)  as  shown  in  figure  1.  The  inside  of  the  shell  .ontains 
a  fluid  of  density  p|  and  sound  velocity  Cj.  In  general,  the  outer  fluid  will  be  different,  and 
is  described  by  parameters  pQ,  C0,  The  elastic' shell  is  a  homogeneous,  isotropic  material  of 
density  p  and  Lame  constants  X,  p.  The  longitudinal  and  transverse  elastic  wave  velocities 
are 

Cl“  =  and  C  j  “  =  ^  ,  respectively. 
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By  symmetry  of  the  incoming  wave  boundary  condition,  we  need  to  consider  only 
the  axisymmetric  problem,  i.e.,  no  variables  will  depend  upon  the  $  coordinate  of  our 
spherical  coordinate  system.  We  will  further  assume  that  there  is  no  motion  of  the  elastic 
shell  in  the  <p  direction 

u*  =  0  .  (29) 

Components  of  vectors  and  tensors  in  this  chapter  refer  to  the  physical  components 
relative  to  the  natural  basis  associated  with  the  spherical  coordinate  system. 

The  displacement  vector  in  the  shell  may  be  represented  in  terms  of  potentials  as 

u  =  70  +  V  X  ij/  .  (30) 

Next,  we  will  demonstrate  that  one  component  of  the  vector  potential  \Ji  is  sufficient  to 

construct  a  solution  to  the  axisymmetric  sphere  problem.  Since  the  potential  \J/  is  not 

unique,  we  can  choose  many  possible  potentials  that  will  give  a  solution  to  the  elastic  equa-  j 

tion.  Recall  that  from  Theorem  A  in  the  Summary  of  linear  Elastodynamics,  there  exists  j 

a  vector  v  such  that  J 

u  =  V2v  .  (31) 

We  can  see  that  a  possible  vector  v  that  will  satisfy  the  symmetry  assumption  u^  =  0  is  I 

chosen  by  letting  v^  =  0.  That  is,  writing  out  the  Laplacian  in  spherical  coordinates  (refer-  j 

ence  5.  p.  1 16), 

u0  =  (V2^  =  , 

r  r-sin“0 

we  can  see  that  =  0  implies  u^  =  0. 

Now,  referring  to  the  proof  of  Theorem  B,  equation  13  says 

\p  =  -  curl  v  . 

In  spherical  coordinates,  with  =  0. 

**0er  +  0e(+^[fr(rvtf)-^J  .  (32) 

Therefore,  a  sufficient  choice  for  the  vector  potential  consistent  with  the  symmetry  of  the 
problem  is  i|/r  =  \pg  =  0,  retaining  only  i//^.  So 

^  =  1^0  £0  • 

Equations  22  and  23  may  be  expanded  in 
(r2<N-)’r  +  sTiTo  (sinO(M’0  +  kL2r2^ 

(f2^, r)  r  +iib(sin0  ’W  o  +  (kT2 


spherical  coordinates  to  give 
=  0 


r^-V) 

cin-fl  / 


^0  =  0. 


(33) 
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Solutions  to  equation  33  may  be  obtained  using  separation  of  van  ’es>  The  radial  equation 
lias  a  solution  in  terms  of  spherical  Bessel  functions  j£  (*),  ng  (*),  and  the  angular  equation 
has  a  solution  in  terms  of  associated  Legendre  polynomials  of  zero  and  first  order,  Pg  (♦)  and 
Pg*  (•),  respectively.  * 


*  =  £g  Pg  (cos0)  Ug  (kLr)  Ag  +  ng  (kLr)  Bj]  = 
£=0 

t 

oo 

^  £g  PfV  (cosfl)  jjg  (kyr)  c£  +  n£  (kyr)  Dgl 


(34) 


where 
and  £g 


arbitrary  constant*  to  be  determined  from  the  boundary  conditions 
included  for  later  convenience.  In  addition,  for  ease  of  writing  let 


Ag1,  Bo,  Cg,  Dg  are  arbitrary  constant*  to  be  d< 
=  (-i)*  (2fi  +  I )  is  included  for  later  conveniei 

Qg  (k[ *)■  ~  £g  Ug  (kgr)  Ag  +  ng  (k|T>  Bg] 
Rg  (kjr)  =  £g  Ijg  (kjr)  Cg  +  ng  (kyr)  DgJ 


Then, 


uu 

0  =  ^  Pg  (e  jsd)  Qg  <kLr) 


*4  =  £  Pg*  (cos0)  Rg  (kyr)  ., 

C=0 

Before  going  any  further,  we  need  some  useful  relations  for  derivatives  of  Legendre 
Polynomials  and  spherical  Bessel  Functions. 

Pg  g  (cost?)  «  Pg3  (cos0)  (37 

Pg  qq  +  cot  0  Pg  q  •■=  -C  (C  +  1 )  1'g  (38 

We  follow  the  sign  convention  of  Magnus  and  Oberhettinger  (reference  6)  for  Legendre 
polynomials..  The  second  relation  follows  directly  from  Legendre’s  differential  equation. 

jgr(kr)  =  kj'p(kr)  (39 

r\rr<kr>  +  2|ifi,r*kr)  +  l*^2"  C(R-M))jg  (kr)  =  0  (40 

The  first  equation  is  a  standard  identity  and  the  second  is  Bessel’s  equation. 
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Writing  out  equation  30  in  spherical  coordinates,  we  find  that 


..  _  ,  ,  1  ,  ,  cotfl  . 


u0  =7^‘  1*4 


(41) 


Using  the  solutions  shown  in  equation  34.  we  find  that 

Po 

P£  + 

e=o 

Using  equations  37  and  38,  we  obtain  the  following  expression: 


"r  “  2  +  ^7%  ♦  ^  V  - 


oo 

“r  '  2  Pt'c“»)[0tT<kLrt- Rflk^)]  . 


e=o 


For  the  radial  functions  let  the  prime  indicate  differentiation  with  respect  to  the  entire 
argument 


HT(kr)  ii  (X)  =  kil'(x) 


ur  =  y  PS(cos0)[kLCrC(kLr)-Mi^R£(kTr)]  . 

e=o 

Similarly, 

OO 

u0  =  J  Pg1  (cos(l)[Qt(rkLr)-  kT  R£  (k-rr)  ~  ~  Rg  (kjf)] 


(42) 


(43) 


£=0 


Expanding  the  stress-displacement  relation,  equation  24,  in  spherical  coordinates,  we  find  that 
T"  =  (In  +  X)  ur  r  +  x[y  ur  +  j  ue  g  +  ^u0]  (44) 

tto  ~  r i  ,  -u0i 
T  -  47  ur,0  +  u0,r  ~  J  * 


(45) 


Substituting  equations  42  and  43  into  equation  44  and  employing  some  algebra,  one  obtains 
the  following  expression: 


Tn  =  y  Pp  (cosfl)  fki  2  ( 2fi Qp  -  XQp)  -  2/jg(g+>  -  (kTRp  -  — )T 


(46) 


Similarly,  equation  45  may  be  written 


T**  =  ji  ^  Pg!  (cos®) 
8=0 


In  the  outer  fluid  media  the  incident  {dace  wave  may  be  written  as  (reference  5. 
p.  1466) 


<W  =  J  * 

e=o 

The  scattered  wave  potential  is  a  solution  to  equation  27  for  outgoing  waves: 

OO 

ds  =  ^  Fj°  Pfi  (cos®)  hg  (k^) , 

8=0 


where 


(‘ 


(' 


hg(*>  =  jg(')  +  ing(-)  .  (± 

The  total  outside  potential  is 

OO 

$0  -  djnc  +  ds  *  ^  $£  Pj  (cos®)  [jg  (k0r)  +  Fg°  hg  (kcr)]  .  <f 

8=0 

Differentiating,  we  obtain  the  outside  displacement: 

uor  =  IT  =  S  P5 (cos0)  tko  j'c  (kor>  +  F8°  ko  h8  <kor)]  -  (i 

8=0 

For  the  stress,  equation  28  gives 

OO 

T£  =  -pQ0i-  J  SgPgtcos®)  [jgtk0r)  +  Fg°hc(k0r)]  .  (! 

8=0 

With  the  requirement  that  the  inside  potential  be  finite  at  the  origin,  we  can  write 

OO 

^  Fg*  Pg(cosO)jg(k|r)  .  (f 


The  displacement  is 


U|r  =  ^  (cosd)  k|  jj  (kjr) .  (55) 

6=0 

and  the  normal  stress 

OO 

T|rr  =  -  pj«-  ^  ^FelPc(c°s0)je(k,r)  (56) 

6=0 

At  the  boundaries  of  the  shell,  we  require  the  normal  displacement  and  stress  to  be 
continuous  and  the  tangential  stress  to  vanish.  Therefore  we  must  satisfy  the  following 
conditions:  At  the  outer  radius,  r  =  a. 


Vx  (r=a)  =  Tan  (r=a) 

(57) 

T^(r=a)  =  0 

(58) 

ur  (r=a)  =  u0r  (r=a) ; 

(59) 

at  the  inner  radius,  r  =  b. 

Trr(r=b)  =  T^lr-b) 

(60) 

(r=b)  =  f. 

(61) 

ur  (r=b)  =  ujr  (r=b) . 

(62) 

For  algebraic  convenienc*.  tel 

x  =  koa  XL  =  kLa  XT  =  kTa 

y  =  kjb  y l  =  kLb  Vj  =  kjb  . 

(63) 

From  equations  46.  53.  and  57  we  obtain  the  following  expression  for  normal  stress 
continuity  at  the  outer  radius,  where  r  =  a: 


-  p0u}~  ^  fyP((cos0)[ij  (x)+  Fg°h2(x)]  = 
6=0 


^  Pg(cosO) 
6=0 


kL"  [-P  0  6<XL^  *  ^Q6^xl0 


2p6(6+l ) 
a 


[■ 


kjR  gUy) 


R((xj)| 


-  m 1UAU, 


Equating  tenn  by  term,  we  find  that 

-  P®"2^  D«  <*> +  f£°  *»£  (*)]  =  kL2  [2nQg  (xL)  -XQ£  (xL)] 

Dividing  through  by  pur  and  using  the  relations 

P  _  P<«>~  .  2  _  ^  . 

?-^  kL 

-7ii[jiW,f(”hiW]  =x4-  f^Qe <*l>  x<?e <*l>] 

28(g-H) 

Xy-  UjRjtxj)*  Rj(x-j-)1  . 

Using  the  definition  cf  Qg  and  Rg,  ghren  by  equation  35.  we  obtain  the  following 
expression: 

[-^hg(x)]  F£0  +  j^;  Upjg  (xL)-  Xjfi  (xL)I  Ag 

+x7v^,J?  *xl^~  ^nc^xL^  ®r  =-^r  IxTjg  (xT)  -  jg  (xth  Cg 

xT- 

_  j-ynj(xj)- ng(xj)J  =  -— ^jg(x) 

Xj“  P 

£>efine  the  following  matrix  coefficients 

S11  _  p  h£(x> 

Sj->  =  0 

513  =  x^rI2pjg(xL)-Xjg(xL)l 

514  =rr>|2/inC(x  L>'XnC<xL)1 

„  _  -26(6+1).  . 

s15 - lxTjg(xT)- jc(xT)J 

XT“ 


(64) 


^16  — {*t) -  "j**!)! 

XT 

Po 

h  =  -  -^-i«  <*> 

Continuity  of  normal  stress  at  the  inner  boundary,  where  r  =  b,  is  similar;  this  time 
we  use  equation  56  for  the  fluid  stress.  By  inspection  we  obtain 

~-jg (y)  f£!  +  iy0 ~  (y !_>* 

l2un'i (yL>-  Xne<yL>l  *t~  lyTJg(yT)-  jg(yT)l  c£ 

yT- 

■  -i  ^  ly-j-ng  (y-p)-  ng(yj)]  Dg  ■=  0  .  (65) 

yT~ 

Define  the  matrix  coefficients 

551  =  0 

Pi 

552 

553  =  X^T  <>L> '  Xi£ 

554  =  xkh  l2i,ne  (yL>  ■  Xn8  iyL>l 
s55  =  -~6  lyjj'g  (yT)  -  jg  (yT)l 

yT- 

%6  =  "  lyj  n«  (>T>  -  nc(>’T^  - 

yr 

From  equations  47  and  58  we  derive  the  following  expression  for  the  tangential 
stress  at  the  outer  radius,  where  r  =  a: 

a  [kLQ£<xL,"~^r1']  "  [kT2Rg(xT,-5*^L>(2-8(«+i)]  =0. 

Using  the  expressions  for  Qg  and  Rg.  we  find 
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21XL*<XL>"  A^  +  2  Ixl^Ul)-  n£<xL)]  B* 

-  *f  'j{  ^XT^ +  I  )J(  Uyjj  Cf 

-  xx2  ri£  (xT)  +  (£+2)  <£- 1 )  n£  (xy)]  D£  =  0 . 

The  matrix  elements  are 

S31  =  s32  =  0 

533  =  2  {xLj£(xL)- j£(xL)l 

534  =  2lxLn£<xL)*  "L^V1 

535  =  ~  (xj)  +  (£+2)  (£- 1 ) j£  (x  j-) J 

536  =  ”  [xT_n£  (xT)  +  (£+2)(C  I)nc(xT)]  . 
Similarly,  at  tne  inner  boundary  r  =  b.  and 

2  fyj/  a£  +  2  i>  Ln£  n£  <>L^  b£ 

-  Jyx~j{  (yx) +  (£+2;  <£- 1 ) j£  tyx)J  c£ 

-  [yx2 n£  (yx)  +  (£+2 )(£-!»  n£  (yx)  J  D£  «  0  . 

The  matrix  coefficients  are 

S4 |  =  S42  =  0 

s43  =  2  •yue^'L*'  h^'O^ 

S44  =  :  |yLn'c(yL*-  "C(yL)l 

s45  =  -  ryX2j{  (yx)  +  (£+2)(£-l)je(y  r)J 

S46  =  ‘  \yj2n'k  (>T>  +  <5+2)  (£-  I )  nc  (yT)  J  . 


Using  equations  42, 52  and  59  we  obtain  tbe  following  expression  for  the  continuity 
of  normal  displacement:  •  - 

*£  Cko»£  (x>  +  ^o1^  <*>]  =  W  <*L>  "  *£  <*T> 

Multiplying  through  by  a,  we  obtain 

£g  [xjG'  (x)  +  Fg°xh'  (x)]  =  xLQf'  (xL)-  £  (£+1)  Rg (xT) . 

Expanding.  we  obtain 

-  xhg'  (x)  Fg°  +  xLjg'  (xL)  Ag  +  xng'  (xL)  Bg 

-  *  (5+1 )  Ijg  (xT)  Cg  +  ng  (xT)  Dgl  =  xjg'  (x> ,  (68) 

Define  the  matrix  elements 
S-*j  =  -  xhg'(x) 

S-»i  ~  0 

S23  =  XL»C  (xL) 
s24  =  xLn5<xL) 

^25  =  ”  5(£+l)jg(xj) 

S26  =  '  5  (5+1 )  ng  (xj) 

^2  =  xJf*  (x)  - 

Similarly,  at  the  inner  boundary  r  =  b.  and 

-  yjg'  (>•)  Fg1  +  yLjg*  (yL)  Ag  +  yLng'  <yLi  Bg 

-  2  (C+l )  |jg’ (yj)  Cg  +  ng  (y-j-)  Dgl  =  0.  (69) 

Define  the  matrix  coefficients 
0 

-  yjg'»y) 

>Lne'(>T) 

-  t  (5+1 )  jg  ;>T) 


•*61 
$62  = 
S63  = 
^64 

S(>5 
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S66  =  -*(«+l)nkryT). 


From  the  boundary  conditions,  we  troy  write  the  following  matrix  relation  for  the 
unknown  constants  A£,  Bg,  Cj,  Dg.  Fg°.  Fg*.  Few  C  >  1 , 


S!1 

0 

SI3 

SI4 

SI5 

s,; 

V 

V 

S21 

0 

% 

*24 

S25 

s26 

F«* 

h 

0 

0 

S33 

S34 

S35 

S36 

A£ 

0 

0 

0 

S43 

S44 

S45 

S46 

b£ 

0 

0 

%2 

S5  3 

S54 

S55 

S56 

0 

0 

S62 

S63 

S64 

*65 

S66_ 

0 

For  fi  =  0,  we  have  a  special  case  since  the  associated  Legendre  polynomials  P£  *  (cosfl) 
which  appear  in  the  tangential  stress  must  vanish.  Therefore,  the  matrix  becomes 


SI1  0  SI3  S14 

F  ° 
ro 

P\ 

S-)  j  0  Sij  Si^ 

F  * 
ro 

Pi 

0  S52  S53  S54 

Ao 

0 

_°  s62  S63  S64__ 

oW 

1 _ 

0 

(71) 


If  we  allow  the  inner  and  outer  fluid  media  to  be  the  same  and  rearrange  the  matrix, 
we  derive  Goodman's  matrix  (reference  4)  except  for  sign  convention  differences.  The 
coefficients  Ag  and  B$>  differ  by  a  minus  sign  from  those  of  Goodman. 


NUMERICAL  IMPLEMENTATION 


In  this  section  the  details  appropriate  to  computer  implementation  of  the  previously 
derived  solution  are  considered.  In  particular,  the  matrix  is  further  manipulated  and  a  meth¬ 
od  of  solution  which  takes  advantage  of  some  peculiarities  in  the  matrix  is  discussed. 

The  standard  recurrence  relations  for  spherical  Bessel  functions  will  be  used  and  .ire 
repeated  here  for  convenience. 

j'(x)  =  xin(x)-jn+l  <x> 


jn"(X)  =  ']  jn(x)  +  xjn+l  <x)- 
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The  matrix  elements  Sg,  after  some  algebra,  may  be  written  in  the  following  form: 

Sll  =’7h*<*) 

S,3  =  [4xL-*t+l  (xL)  +  (2fi(8-l)-  xT2)  jg(xL)] 

S14  =  ^2  [4xL**«  <xL)  +  (2*(«-l)- xt2)  n£(xL)] 

Sj5  =  -2*<*+1)  I(C-l)jc(xT)-  xTj£+,  (xT)J 
e  _  -2£f£+I) _ 

16  ^2  f(8-I>n£(xT>-  xTn£+l  <xT)l 

S21  =  -fih£(x)  +  xh£+I  (x) 

523  =  {ic<XL>-  xLi«+l  <XL> 

524  =  knfi(xL>"  xLn£+l  (xL> 
s25  ~  -®  (f+l)jjtxT) 

S26  =  “*  <k+Dn£(xT) 

1  (72) 

s33  =  2  ICC- 1  >  jg  (xL)-  xLj£+I  (xL)l 

534  =  2l«!-l)n£(xL)-  xLn£+1  (xL)l 

535  ~  [XT~  “  3  (2~-  1)J  j£  (Xj)  -  2xjj£+j  (xT) 

%6  ~  [*r'  “  2(C“-1)J  ng(xj)-  2xjn£^j  (xj) 

543  =  2  UC- 1 ) j£  (yL)  -  yLjg+j  (yL)] 

544  =  2K«-l)n£(yL)-yLr.jH  (yL>) 

545  =  [yT2  -  2  (fi2-  1 )]  jfi  (yT)  -  2yTj£+,  (yT) 

546  “  [yT  “  2(C--])J  n£(yT)-  2yjn£+I  (yy) 

S52  =  7^  (y> 
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jy  l .  • 


c.  -  1 

s53 - 2 

y-r 

[4yLje+i  (yL> +  (m(«-D-5T2’ 

1  ie  (vL)] 

S54  =“4 
yT 

4yLne+i  (yL)+  (2E (E- 1 ) -  yT2) 

1  "E  (EL>J 

S55  =  — [(C-  l)jg(yj)-  yjjg+j  (yj)l 
yj 

%6  =  -  ^  ^  {(fi- 1 )  ng  (yj)  -  yjng+j  (y^)] 

VT 


s62  =  -%(y>  +  yjg+i  (y) 
s63 =  %  (y^  ~  yuc+i  W 

s64  =  fing(yL)-  yLng+1  (yL) 


(72  contd) 


S65  =  -C(C+l)jg(yT) 
=  -C  (C+I)  ng  (yj) 


h  =  yMx> 

@2  ~  cJg(x)“  xJg+j  00  • 

Notice  fhat  only  two  matrix  elements  are  complex..  S  ]  ]  and  j .  Rather  than  use  a 
complex  matri\  solution  routine,  a  solution  algorithm  is  discussed  which  will  allow  one  to 
construct  the  solution  of  the  con:ole  <  matrix  by  solving  a  real  matrix  with  two  different 
right-hand  vectors. 

Divide  the  first  column  of  the  Sy  matrix,  equation  70,  by  S2 1  - 


S1  l/S2I 

0 

S13 

S14 

S1 5 

S16 

VS21 

~P\ 

1 

0 

s23 

S?4 

S25 

*26 

Ffi1 

h 

0 

0 

S33 

S34 

S35 

S36 

Ae 

0 

0 

0 

S43 

s44 

S45 

S46 

Bg 

0 

0 

S52 

S53 

S54 

S55 

S56 

0 

0 

S62 

S63 

S(i4 

S65 

s66 

Dg 

0 

(73) 


Written  in  vector  form, 

Sz  =  j3. 

Expanding  in  terms  of  real  and  imaginary  parts,  we  find  that 
(SR  +  iSj)  (zR  +  izj)  =  0  : 

or 

SR~R "  SI^I  =  a 

SR  ~I  +  SI  i R  =  0  • 

Let 

“r  ■ Re  (St)  •“! =  lm  (Si)  • 

Then,  the  real  and  imaginary  parts  of  the  matrix  S  are 


“R 

0 

S13 

S14 

S15 

S16 

1 

0 

S23 

S24 

S25 

S26 

0 

0 

S33 

S34 

S35 

S36 

0 

0 

S43 

S44 

S45 

S46 

0 

S52 

S53 

S54 

S55 

S56 

0 

S62 

S6? 

S64 

S65 

^66 

The  imaginary  part  has  only  one  element: 

SI  =  “i £l  ©3|  ■ 

Using  this  fact,  equation  76  may  be  written 

SRil  =  '  S1  ^  =  "  al  ZR1  £|  • 

Suppose  we  solve  the  following  equation  for  a  vector  x. 

SR  £  =  ~1  ' 

Then,  it  is  clear  that  zj  is  some  scalar  multiple  of  x,  so  we  write 
l\  =  . 


where  the'  constant  7  is  to  be  determined. 


where 


“R 

0 

S!3 

S14 

S15 

S16 

¥ 

0 

0 

$23 

S24 

$25 

$26 

Ftl 

0 

0 

S33 

S34 

S35 

S36 

a£ 

0 

0 

S43 

S44 

S4S 

S46 

Bg 

0 

S52 

S53 

S54 

S55 

S56 

C£ 

G 

S62 

S63 

S64 

S.65 

S66_ 

^23  = 

S13" 

S23  aR 

^25  = 

S1S" 

S25aR 

s24  = 

S14" 

S24  aR 

B26  = 

$16  ~ 

S26  “R 

*1 

h 

o 

0 

0 

0 


¥  -  ¥/sn 

02  =  (?]  -  02  aR  • 

Eliminate  the  first  row  and  column  of  equation  89. 


0  $23  $24  S25  S26 

"V 

02 

0  S33  S34  $35  $36 

AC 

0 

o  s43  s44  s45  s46 

Bfi 

= 

0 

S52  S53  S54  S55  S56 

ce 

0 

_S62  S63  S64  S65  $66_ 

3- 

0 

and 


aR¥  +  s13AC  +  S14BC  +  Sl5C5  +  S16DE  =  Pi  • 


(89) 


(90) 


(91) 


(92) 


(93) 


row. 


Similarly,  multiply  the  sixth  row  of  equation  92  by  -  S52/S62  and  add  to  the  fifth 
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0  S23  S24  S25  ^26 

Fe* 

H 

0  S33  S34  %5  S36 

* 

0 

o  s43  S44  S45  s46 

Bg 

= 

0 

S52  S53  S54  S55  S56 

Cg 

0 

_  0  S63  S64  s65  S66_ 

3- 

0 

Then,  eliminate  the  first  column  and  fifth  row. 


% 

B24 

B25 

S26 

a£ 

h 

S33 

S34 

S35 

S36 

B£ 

* 

0 

S43 

S44 

S45 

S46 

0 

1 

J*3 

%4 

hs 

B66_ 

_°£. 

0 

?52FCf 

+  S53A£ 

+  S54 

B£  +  S55C£  +  S56 

Dg  = 

(95) 


(96) 


The  real  matrix  solution  is  therefore  reduced  to  the  solution  of  a  4X4  matrix  and 
the  appropriate  multiplications  to  find  Fg^and  Fg°  from  equations  96  and  93,  respectively. 

As  we  saw  before,  the  fi=0  math;:  is  a  special  case  of  equation  71 .  We  can  follow 
the  same  procedure  outlined  above  to  solve  the  full  matrix.  The  real  matrix  becomes 


"R 

0 

S13 

S14 

1 

0 

S23 

S24 

1 

0 

S52 

S53 

S54 

0 

S62 

S63 

S64 

(97) 


We  may  also  follow  a  similar  Gaussian  elimination  procedure  to  take  advantage  of 
the  zeroes.  Following  the  same  steps  as  before  we  must  solve  a  2X2  matrix. 


r_  _  1 

r  1 

1  1 

s23 

S24 

> 

0 

11 

h 

_  B63 

%4_ 

I - 

CC 

O 

1 _ 

0 
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*V»^1 


-71535255 


- r-TT  '”'  * 


where 


S23  ~  Sl3_aRS23*S63  "  S53"  %3 


S24  "  SI4"  aRS24  S64  "  S54"  S* 


°RFo°  +  S13Ao  +  SI4Bo  ~  h 


(100) 


S52Fo1  +  S53ao  +  S54Bo  =  0 


(101) 


In  conclusion,  the  lime-saving  methods  described  in  this  section  have  been  incorpo¬ 
rated  into  the  spherical  shell  computer  program.  Comparison  of  a  previous  version  of  the 
program  using  a  standard  LU  decomposition  of  the  original  6X6  matrix  and  the  modified 
version  shows  a  60-percent  savings  of  the  actual  inversion  time,  resulting  in  an  overall 
computer  program  time  savings  of  35  percent. 

LIMITING  CASES 


We  will  now  consider  several  limiting  cases  of  the  elastic  shell  solution.  In  the  limit 
of  vanishing  shell  thickness,  the  solution  is  shown  to  reduce  to  scattering  from  two  dissimi¬ 
lar  fluid  media.  In  particular,  when  the  wavelength  is  much  greater  than  the  size  of  the  fluid 
sphere,  we  obtain  a  solution  applicable  to  such  Ixvdies  as  air  bubbles.  If  the  shell  thickness 
is  non-zero  but  sufficiently  thin,  a  first  order  approximation  to  the  solution  is  obt ?*!  <  <( 
using  a  Taylor  series  expansion  of  the  matrix  elements.  This  solution  is  useful  to  avotd 
linear  dependence  problems  in  the  Sjj  matrix  for  sufficiently  thin  shells. 

Suppose  we  let  the  thickness  of  the  shell  l>e  zero,  a  =  b..  T  hen.  by  inspection  of  the 
matrix  elements,  equation  72. 


Slj  "  S5j 
S2j  _  **6j 
S3j  “  S4j 


j  =  3. 4.  5. 6 


(102) 


Subtracting  the  fifth  row  of  the  shell  matrix,  equation  70.  from  the  first  row  and  the  sixth 
row  from  the  second,  we  obtain  a  relation  wl’  ch  determines  the  coefficients  for  the  inner 
and  outer  fluid  potentials: 


~S52~j  h°l  _  [01 


!l  "°62 


-SAs  I 


kl- 


(103) 
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%f,.- .  xV\S*  *4  •At*** W* k* 


Ff°l  r^w] 

-  xh |>(x)  Ff*  J  li^W  J 


(10*) 


-  *t  -  Co 

•bone  n  =  T~  ~  ~F~ 


The  inverse  of  the  matrix  is  easiy  obtained: 

[Ff°n  |  r+ipjf<nx)  +/»jji(*ix)j  r-Poitu, 

Ff(l,J  [  **4 «  P0h f(x)j  [  x£(x) 

•fi-eie  A  =  *  p0ip&£ (x> jj (ijx) -  pjxhg(x)  jj (ijx)  . 

The  coefficient  for  the  outer  polenta!  is  easily  obtained: 

0  _  -  Pnq>f  (x)ip(ipO-»-Pfjf  Cx)if  |nx) 
rfc  +  PoW <*)  J£  (*!*)  -  Pf1^  'x)i  <r/*)  * 


(105) 


where 


o  =  [jg  (x)  +  iT  gjt  (x)  1 

8  |h£(x)  +  ir2hg(x)J  ’ 

»P0r/  fjg  fax)1 

c  ~  P|  Lifi(7»x>J 


(106) 


This  result  is  identical  with  Morse  and  Ingard  (reference  8)  except  for  sign  convention. 

When  the  wavelength  is  much  greater  than  the  size  of  the  sphere,  we  may  use  a 
small-argument  limiting  form  of  spherical  Bessel  Functions  (reference  9) 


j{(x) 


ng(x) 


In  particular. 


1.3 .» (2C+1 ) 

- 1.1.3 -(2g-l] 


x  -*•  0 


2* 


"'SS!^H85^53^Kr" 
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Droppiu  higher  order  terms  in  denominator. «  find 


The  hither  order  coefficients  are  proportional  to  x~®+*  and  are  therefore  negligible. 

At  the  surface  of  the  sphere  we  can  obtain  an  expression  for  the  total  potential 
using  the  first  two  terms  of  equation  5 1 . 
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Notice  that  for  a  light.  compressible  sphere  like  an  air  bubble,  where  P|  «  pc.  the  linear 
term  in  x  is  neslieiblc.  Therefore. 
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Tim  neurit  ^rces  whir  that  obtained  for  acoustic  scattering  from  air  babbles  m  »a!cr  (refer- 
eacesSO.ll). 

Whea  the  thkkncsofthe  shd  b  ssssaH.  bat  aoe-mo.  a  fnt  order  2pprouiut>oo  to 
Ok  tfcdl  matrix  may  be  ohedaed  uaqg  a Taylor  series  expanrioo  of  the  matrix  dements. 

Tk  procedure  does  not  provide  much  wlyticil  srinpiifkation  of  the  matrix,  but  avoids 
numerical  diffkuiticf  in  solution.  When  the  did  is  linn,  two  rows  of  the  shell  matrix  are 
approximately  equal,  creating  a  thready  dependent  matrix. 

la  order  to  corspkmcnt  the  numerical  solution  method  described  in  the  previous 
section  for  thin  shells,  we  will  wort  with  the  reduced,  real  matrix,  equation  95. 

$23  ^24  ^25  *26  *t  h. 

%3  ^4  %6  S36  ®t  0 

=  (112) 
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where  the  barred  matrix  elements  are  defined  by  equations  90  and  94. 

Notice  tint  for  a  thin  shell,  the  second  row  is  approximately  equal  to  the  third  row  of 
the  matrix  eqiution  111.  Expanding  the  elements  in  the  third  row  in  a  Taylor  s  series  about 
«=i  wc  obtain,  to  first  order  in  the  thickness  5.  these  expressions  for  the  matrix  elements: 

S43<YL*  =  S33(xL>~  kLfiS33(xL> 

S44(YL*  ~  S34(XL,_ 

(113) 

s45  (YT)  “  S35  (XT) '  kT5S35  (xT) 

S46  (YT)  =  S36  (XTJ "  kTSS36  (xT> 

Subtracting  row  3  from  row'  4.  we  find  that  the  matrix  equation  1 1 2  becomes 
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The  matrix  is  valid  to  first  order  in  the  thickness  6  and  eliminates  the  linear  dependence 
problem  when  the  shell  is  sufficiently  thin. 
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In  tins  section  typical  computed  results  are  presented  for  scattering  from  an  xr-flkd, 
itaMMU  spherical  shell.  The  results  shuva  in  figures  2  and  3  were  obtained  using  the  fol¬ 
lowing  input  data.  MKS  units  uc  used  throogboul. 

p{  =  1.21  C,  *  343 

l>0  =  1.026X103  C0  =  1.5X103 

p  *  2.7X103  CL  =  6.4I2X103  Cj  =  3J043X103 

i/a  =  .05  r/s  =  I 

The  total  pressure  as  a  function  of  at  the  surface  of  the  elastic  shell  is  shown  for 

angles  from  the  incoming  wave  of  0  and  180  degrees,  respectively.  The  complexity  of  the 
scattered  field  is  due  to  the  vibrational  characteristics  of  the  elastic  sphere. 


Predicted  total  pressure  as  a  function  of  k0a  at  the  surface  of »  spherical  shell,  0  •  1 80  degrees. 
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